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We characterize the heterogeneous character of the dynamics of liquids approaching the glass tran- 
sition through an experimental determination of the number of dynamically correlated molecules 
Neon- as obtained from dynamical susceptibilities. To do so, we have obtained a new set of dielectric 
spectroscopy data for liquid dibutyl-phtalate on a fine and extended temperature and pressure grid, 
and we have used it in conjunction with high-pressure data from the literature. We have been able 
to evaluate the contributions to N CO rr that are due to fluctuations associated with density and with 
temperature separately, thereby improving the estimate of N CO rr- We find that N corr increases along 
the glass transition line, and more generally along any isochronic line, as pressure increases (up to 
1 GPa), a result which is at odds with recent reports and theoretical predictions. 



I. INTRODUCTION 

A tempting explanation of the dramatic slowing 
down of the dynamics observed in liquids approaching 
their glass transition is the existence of a growing un- 
derlying length scale. The relevance of such a length, 
characterizing the cooperative nature of the structural 
relaxation, was already suggested by Adam and Gibbs 
in the sixties [1]. The search for length scales went in 
two directions [2|: on the one hand, studies of non- 
trivial spatial correlations in the structure, associated 
with some hidden order parameter or in the form of 
point-to-set correlations describable through the influ- 
ence of amorphous boundary conditions 1 3 1 (both being 
essentially inaccessible in laboratory experiments); on 
the other hand, based on the mounting evidence in ex- 
periments and simulations of the increasingly heteroge- 
neous character of the dynamics [4-6 |, investigations of 
a "dynamic" length describing the spatial extent of the 
dynamic heterogeneities. However, direct experimental 
measurement of the latter in molecular liquids is hard 
to obtain. In consequence, the corresponding results are 
scarce |4j [8j, moreover providing no information on 
the evolution with temperature. 

A few years ago, on the basis of theoretical arguments, 
a new method for estimating the typical length scale 
of the dynamical heterogeneities and access its temper- 
ature dependence was suggested It relies on the 
introduction and the study of multi-point dynamical 
susceptibilities. No information on spatial correlations 
in the dynamics can indeed be derived from the mere 
consideration of standard 2-time correlation functions, 
F(f) =< SO(0)SO(t) > (with O(f) some observable de- 
pending on the probe), which describe the "average" 
dynamics of a liquid. Higher-order space-time correla- 
tion functions, or their volume integrals that represent 
susceptibilities and may therefore be more easily acces- 
sible, are required to probe the fluctuations around the 
average dynamics, whose manifestation is precisely the 



dynamic heterogeneities |10|. From the multi-point dy- 
namic susceptibilities, one can extract an estimate of the 
"number of dynamically correlated molecules" ||9j [TT1 . 
which is the focus of the present work. 

Specifically, the four-point dynamical susceptibility 
Xi(t) = N < 5C(0,t) 2 >, where N is the number of 
molecules in the system and SC(0, t) = SO(0)SO(t) — 
F(t) characterizes deviations from the average dynam- 
ics, quantifies the overall amount of spatial correlations 
between spontaneous fluctuations taking place between 
times and t. This quantity is still, unfortunately, very 
difficult to directly measure in experiments, especially 
in the case of molecular liquids. Interestingly, how- 
ever, it is related through fluctuation-dissipation rela- 
tions [9, 12J to other dynamical susceptibilities associ- 
ated with fluctuations that are induced by an external 
control parameter x (such as temperature, pressure or, 
for colloids, volume fraction): these are three-point sus- 
ceptibilities, defined as Xx{t) = / which describe 
the response of the two-time correlation function F(t) 
to a change in an external parameter. Such quantities 
are much easier to experimentally access than Xi(t)> as 
the dependence of F(f) on temperature, pressure, or 
any external field, may be measured by several well- 
developed methods such as dielectric or photon corre- 
lation spectroscopy for example [9]. It was moreover 
shown by computer simulation studies [12] that a good 
estimate for Xi(t) is obtained from the Xx{t)'s in the 
viscous liquid regime. (Note that another 3-point dy- 
namical susceptibility has also been directly measured 
by means of an investigation of the nonlinear dielectric 
response 1131.') 

Thanks to these high-order dynamical susceptibili- 
ties, direct experimental evidence of the increase of the 
number of dynamically correlated molecules N corr (T) 
while approaching the glass transition by cooling was 
obtained for several molecular liquids and polymers 
IlIU EMI. The growth of N corr (T) was found to be 
moderate compared to the dramatic change in relax- 
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ation time in the same temperature range, with JV corr 
reaching a value of the order of 10 2 at the glass transition 
temperature To at atmospheric pressure. No clear con- 
clusions could be drawn concerning the influence of the 
chemical nature of the system on the variation of N CO rr- 
In particular, no correlation was observed between the 
growth of N corr (T) and the fragility of the system, i.e. 
the deviation of the temperature dependence of relax- 
ation time from an Arrhenius behavior. 



When cooling a liquid at constant pressure, its den- 
sity increases. It is therefore useful to disentangle the 
effect coming from the densification of the system from 
that due to the decrease of thermal energy. This pro- 
vides clearer insight into the physics of the slowing 
down of relaxation JTHJ. In addition, when consider- 
ing estimates of the four-point dynamical susceptibility 
Xi(t), it may prove important because a better descrip- 
tion of the associated fluctuations in the NPT ensemble 
is provided by summing the separate contributions of 
the density-induced and the temperature-induced fluc- 
tuations. However, except for partial results at atmo- 
spheric pressure ITU , this latter procedure was never 
attempted. Previous pressure studies along several iso- 
bars considered the crudest bound on Xi (t ) only II161I17I . 

In this work, we consider the effect of density and 
temperature on the number of dynamically correlated 
molecules in a fragile glass-forming liquid, dibutylphta- 
late. To do so, we have measured the (linear) dielectric 
susceptibility along both isobars and isotherms, cover- 
ing thermodynamic points on a fine grid in temperature 
and pressure, and combined these results with existing 
high-pressure data [19 20 1. We have calculated the re- 
sponses of the linear dielectric susceptibility to tempera- 
ture along isobars and to density along isotherms to pro- 
vide estimates of the variation of the number of dynami- 
cally correlated molecules. To go beyond this and, as ex- 
plained above, improve the estimate of N corr by consid- 
ering separately the contributions associated with den- 
sity and temperature at any state point, we have used 
the density-temperature scaling already shown to de- 
scribe relaxation data in supercooled liquids in a large 
domain of temperature and density [21-25]. With the 
help of this scaling description, one can more easily 
compute the various dynamic quantities entering in the 
expression of N corr . We have taken special care to min- 
imize and to evaluate the systematic errors that come 
with the calculation of Ncorr- Our main finding is that 
Ncorr significantly increases (by a factor of almost 3) 
along the glass transition line as one increases the pres- 
sure up to 1 GPa (the density changes by 25%). The 
same trend is found for other isochronic lines (i.e. lines 
of equal relaxation time). Contrary to previous reports 
lElHZI and to theoretical predictions [26, 27J, the extent 
of the spatial correlations in the dynamics is therefore 
not uniquely determined by the relaxation time. 
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FIG. 1. Imaginary part of the dielectric spectrum of liquid 
DPB measured at T = 216 K from atmospheric pressure (left- 
most curve) to P = 250 MP a (rightmost curve) obtained by 
both compression and decompression. The average pressure 
step between 50 and 250 MP a is about 10 MP a. 



II. EXPERIMENTAL DETAILS 

A. Sample and experimental setup 

The dibutylphtalate (DBP) sample (99% purity) was 
acquired from Sigma Aldrich and used as purchased. 
The experimental pressure dielectric setup for the di- 
electric measurements was described in a previous work 
1281 . The main characteristic of this setup is that it in- 
sures a completely isotropic compression by the use of 
an external pressure fluid. For the present study, we 
measured a new set of data with very small pressure 
or temperature steps (see Fig. [T] for an illustration). 
This data set is fully compatible with that reported in 
Ref. [28 j but is more complete in order to perform the 
detailed analysis involved in estimating N corr . The re- 
producibility was always checked very carefully, espe- 
cially under pressure where the experiments were car- 
ried out alternatively by compression and decompres- 
sion for comparison. 

The dielectric spectra were fitted by Havriliak- 
Negami expressions 11291 : 



e(w) 



Ae 



[1 + (icor^p ' 



(1) 



with Ae, t, ol and 7 temperature and pressure depen- 
dent adjustable parameters. This provided a convenient 
and faithful parametrization of all data, allowing for 
instance to conveniently extract the relaxation function 
through Fourier transformation and to take numerical 
derivatives with respect to the control parameters via 
the dependence of the adjustable parameters. 
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B. Thermodynamic data 

Information on the thermodynamics is needed for the 
computation of the number of dynamically correlated 
molecules and it must be considered with care in or- 
der to minimize the sources of uncertainty. The PVT 
data was obtained from the literature and included very 
high pressure results [19.. 30, 31 1. (This has been already 
described in Ref. [28 1.) The density and temperature 
dependent isothermal compressibility Kj(T,p) and the 
thermal expansion coefficient a.p(T,p) were estimated 
from this data set. 

Some results for the temperature dependence of the 
heat capacity Cp of DBP are available in the literature 
11321 |33l , but only at atmospheric pressure. We esti- 
mated from thermodynamic arguments and compar- 
isons to other molecular liquids that the pressure vari- 
ation of Cp is negligible for DBP. This is in part due 
to the fact that the absolute value of Cp is very high 
for DBP (about 440 /K _1 moZ _1 at Tg and atmospheric 
pressure [32J). The variation for the constant volume 
heat capacity Cy can in turn be estimated from the re- 
lation Cy = Cp — Tap V/kt and is also found to be 
small (typically less than 10%). The jump of Cp at 
Tg is known at atmospheric pressure [32 1 and is about 
150 JK~ 1 mol~ 1 . Unfortunately, no data are available for 
the jump of Cy nor for that of Kj at the glass transi- 
tion. (No data either are available for the pressure de- 
pendence of these jumps.) For other molecular glass- 
forming liquids, one observes that the heat capacity Cp 
of the liquid at TV slightly increases with pressure: a few 
% for 3-methylpentane and 1-propanol up to 200 MP a 
1341 , a few % for m-fluoroaniline up to 400 MPa (351, 
and 10 - 12 % for toluene up to 400 MPa |36l- However, 
the Cp of the glass at Tg has been found to also increase 
in 3-methylpentane and 1-propanol where it has been 
measured [34], so that the heat capacity jump ACp ig (P) 
barely changes with pressure in these systems. 



III. NUMBER OF DYNAMICALLY CORRELATED 
MOLECULES N corr ALONG ISOBARS AND ISOTHERMS 

We are interested in evaluating the four-point dynam- 
ical susceptibility X4 an d the associated number of dy- 
namically correlated molecules in the NPT ensemble. 
From considerations about fluctuation-dissipation rela- 
tions [9 , 12J, two expressions may be used. The first one, 
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where cp is the isobaric heat capacity per molecule, re- 
lates X4 FT (t) t° tri e three-point-dynamical susceptibil- 
ity Xr PT (t) = (ttt^) ' W1 th a strictly positive residual 
term describing the fluctuations in the NPH ensemble 
(where H is the enthalpy). The second expression takes 



into account the contributions associated with tempera- 
ture (at constant volume) and density (at constant tem- 
perature) separately: 



xr T (t) 



k B T 2 , 
cv 



NVT (t)] 2 +p 3 k B TK T [ X f T (t)f- 



(3) 



where cy is the isochoric heat capacity per molecule, 

j#vt(o = (W) p ,x? VT v = (^) T '» d *r E is 

the (strictly positive) four-point-dynamical susceptibil- 
ity in the NVE ensemble. For a more detailed discussion 
concerning these different quantities, see Refs. IfTTI and 

na. 

In simulation studies on model glass-formers |12|, 
Xa an d X^ VE were found to be small compared 
to the other terms in the viscous liquid regime at 
low temperature, becoming even negligible at the 
lowest temperatures. This was also indirectly con- 
firmed by a comparison between results obtained from 
(k B T 2 /c P )[xT PT (t)} 2 and from the direct experimen- 
tal measurement of a nonlinear dielectric suscepti- 
bility 1 13 J . It should nonetheless be stressed that, 
since Xi VE < Xi > considering the separate con- 
tributions associated with temperature and density, 
(k B T 2 /c v )[x^ VT (t)} 2 + p 3 k B TK T [xf T {t)] 2 , provides a 
better approximation to the dynamical four-point sus- 
ceptibility than {k B T 2 /c P ) [x^ PT {t)] z . 

In this work we focus on the numbers of dy- 
namically correlated molecules N corr which may be 
defined as the maximum over time of the rele- 
vant time-dependent susceptibilities (or their absolute 
value if negative). This maximum takes place for a 
time of the order of the average relaxation time of 
the system t k . One therefore considers N^ f r T 4 = 
max t [x^ PT (t)}, NZ T T = Vks/cp Tm«x t { |^ PT '(f) I }, 



N NVT T 
corr,T 



Vk^TmaxtHx^m^and N™? 



y/p3k B TK T max t {xP PT {t)}. From Eqs. § and g, N^ 4 
may then be estimated either from 

N NPT A > (N NPT T ) 2 (4) 

or from 



N NPT > ( N NV T ) 2 , (N NPT n2 



(5) 



with the latter expression providing a better approxima- 
tion than the former one. 

We have first evaluated the four-point dynamical sus- 
ceptibility along isobars through the simpler formula 
in Eq. To do so we have computed Xt PT {t) 

from a numerical derivative of the Fourier transform of 
the Havriliak-Negami fits of the experimental dielectric 
data, as explained in detail in Ref. [TTT1 . We display 
the resulting estimate for N~l*Ji versus temperature T 
on different isobars in Fig. [2] (a). The results are simi- 
lar to those obtained at atmospheric pressure in previ- 
ous papers IflTl 1141 Il5l and in Ref. Ifl6ll : the number of 
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ing section. Note finally that the span of most curves 
in Figs. 2 (a) and (b) is quite restricted. In spite of the 
large number of state points studied and the fine T, P 
grid covered, the range of relaxation times that are di- 
rectly accessible experimentally on given isobars and 
isotherms remains limited, especially at low tempera- 
ture and high pressure. 



IV. N corr THROUGH THE SCALING DESCRIPTION OF 
THE RELAXATION TIME 



We have already stressed that considering separately 
the fluctuation contributions due to density and tem- 
perature provides a crisper estimate of the number of 
dynamically correlated molecules N^. However, to 
compute the various terms required [see Eqs. l|3| and 

f], an alternative procedure to the direct differentiation 
the experimental spectra must be employed. We have 
therefore relied on the, by now, well established scaling 
description of the density and temperature dependences 
of the relaxation time. 



A. Density-temperature scaling law for relaxation time 



FIG. 2. -a Temperature dependence of the estimate of N^fJ, 
from Eq. (jij along several isobars, -b Density dependence of 
N^PTp along several isotherms : closed squares T = 216 K, 
closed triangles T = 220 K, open squares T = 227 K, crosses 
T = 295 X (the data at T = 295 K are from Ref. (20)). Both 
quantities have been computed by a numerical differentiation 
of the Havriliak-Negami fits of the dielectric spectra. 

dynamically correlated molecules is found to increase 
when the temperature is decreased along any isobar, in- 
dicating growing spatial correlations in the dynamics as 
one approaches the glass transition. 

We have also considered the contribution coming 
from the density-triggered fluctuations along isotherms, 
i.e. Nwrr;p(p)- To compute this term, we have nu- 
merically differentiated with respect to density the 
Havrilliak-Negami fits (after converting the pressure 
data to density data) and included the necessary ther- 
modynamic factors. The results are shown in Figj2|(b). 

We observe that N^^F„ increases with increasing den- 
sity along any isotherm. The other contribution involv- 

iftg/tr^W = ( ) j * s virtually impossible to get di- 
rectly from the dielectric data with a good accuracy, due 
to the fact that data are not collected along isochores. A 
different method will then be considered in the follow- 



The relaxation times T a obtained from the Havriliak- 
Negami fits of the dielectric spectra [r a is equivalent to 
t in Eq. ([TJ] are plotted versus pressure at different tem- 
peratures in the inset of Fig. p](a). The main panel of 
this figure shows that the same points, combined with 
literature data [19.. 20 1, can be collapsed onto a master 
curve when plotted versus a scaling parameter e(p)/T: 
the function e(p) is obtained from the best collapse of the 
whole data set and has the physical meaning (up to an 
overall constant factor) of a bare or noncooperative acti- 
vation energy characterizing the dynamics at high tem- 
perature Bill . The fact that all the t^'s can be plotted on 
a master curve, 

log(^)=^), (6) 

with T a species-dependent scaling function, has been 
already proposed and verified for a great variety of 
glass-forming liquids and polymers [21-25J. In many 
cases, the function e(p) can be approximated by a power 
law, p x , [22-24J. Often, however, the density range un- 
der study is too limited to allow a discrimination be- 
tween different functional forms for e(p). Interestingly, 
in the case of DBP, data are available up to 1 GPa, and 
it was found in Ref. |28| that a simple power law does 
not provide a good collapse of the data over the whole 
density range and should be replaced by a more general 
function, e.g. a polynomial function. This is illustrated 
in Fig|3](b). In such a case, the exponent x can be gener- 
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alized to a density dependent parameter 



x(p) 



dlog(p) 



(7) 



which varies between 1.5 and 4 for DBP over the range 
of density covered. As discussed elsewhere[37| (and re- 
cently confirmed in the case of "strongly correlating" 
liquids |38|), there are no physical reasons for why the 
density-temperature scaling law of the relaxation time 
should always be expressible in terms of a simple power 
law. In this respect, the study presented in this paper is 
more general than previous ones Ifl6l IT7I that have only 
considered liquids and ranges of pressure for which the 
exponent x can roughly be taken as constant. We will see 
that this feature has some important consequences on 
the variation of the number of dynamically correlated 
molecules. 



B. Density and temperature contributions to N 



NPT 
corr,4 



The density-temperature scaling of the relaxation 
time provides a continuous parametrization of the lat- 
ter over the whole domain of density and tempera- 
ture and more specifically allows us to compute deriva- 
tives with respect to the temperature at constant den- 
sity an otherwise horrendous task experimentally. To 
further facilitate the computation of the various con- 
tributions to N^?Ja, we have simplified the form of 
the relaxation function F(t) by fitting the result of the 
Havriliak-Negami parametrization to a stretched expo- 
nential, F(f) = exp[— (t/T a )P], with f> a state dependent 
parameter. This enables us to rewrite N^J r T T , N^ ^ p , and 
N NPT 



corr,T 



as 
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NPT 
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FIG. 3. -a Logarithm (base 10) of the (dielectric) relaxation 
time T a (in units of tq = lsec) of DBP measured at differ- 
ent pressures and temperatures versus the scaling parameter 
e(p)/T. The latter is adjusted to obtain the best collapse of 
all data points. The blue symbols denote our new set of data, 
the black ones our previous data from Ref. 1251 , and the vio- 
let ones are for literature data from Refs. |20|and[30]|. Inset: 
All our data points plotted versus pressure, -b Bare activation 
energy e(p) obtained from the collapse of relaxation time data 
(dotted line). For DPB the function cannot be simply described 
by a power-law dependence p x : the full line represents the best 
fit with a constant x = 2.5. 



and 



N. 



NPT 
corr,T 



fie 



1 + x{p)Ta.p 
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where dp is the thermal expansion coefficient, f> is the 
stretching parameter, is the inverse of Euler's math- 
ematical constant, which should not be confused with 
the activation energy e(p), x(p) is defined in eq. and 
J-' is the first derivative of the scaling function in Eq. 
|[6). To derive the above equations, we have used the fact 
that the maximum over time of the various susceptibili- 
ties calculated from a stretched exponential description 
occurs for t — r a . 



Finally, N^ r \ is estimated from Eq. H as 



N, 



NPT 



corrA 



(11) 



(10) where N^ r T T is given in Eq. &8\ and R is defined as IfTTI 



_R = pTKjCy x(p) z 



(12) 



R is therefore a measure of the relative importance of the 
density-induced contribution to the number of dynami- 
cally correlated molecules at constant pressure. 

Before discussing the results obtained from Eq. (TT) , it 
is worth assessing various sources of uncertainty in the 



calculation. We have thus compared the values of 
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and of N^^j computed from Eqs. (9| and 1 10 \, respec- 
tively, to their counterparts obtaineclthrougna numeri- 
cal derivative of the Havriliak-Negami parametrization 
of the dielectric data (see previous section). We find that 
the difference between the two sets are at most of the or- 
der of 10 % [which amounts to at mos t 20 % in the square 
values entering in N™^: see Eqs. 1 4 5 1], with no sys- 



tematic trend with density or temperature. 

Additional errors come with the account of the ther- 
modynamic input. First, there can be a change of an 
overall factor for the absolute value of the N CO rr's accord- 
ing to whether one considers the total fluctuations of 
energy, enthalpy or density (as rigorously enters in the 
lower bound for N™£fy or only that part of the fluctu- 
ations which is associated with the structural relaxation 
and can be ascribed to the values in excess to the glass 
ones (as physically motivated) ITTTI . For instance, as de- 
scribed in section II-B, ACp is about 1/3 of Cp for DBP 
at the glass transition, which, when replacing Cp by ACp 
in the estimate of the number of dynamically correlated 
molecules, introduces an overall multiplying factor of 3. 
However, we are not interested in absolute magnitudes 
but in relative variations with T, p, or P. Concerning the 
latter, the situation is much more favorable, as one can 
estimate that the relative change in the various thermo- 
dynamic factors is never more than 10 % over the whole 
range of state points under study. To summarize this 
discussion, we conservatively conclude that a variation 
of the estimate of N%£F A that is beyond 20 - 30 % is phys- 
ically significant. Smaller ones need to be interpreted 
with caution. 

Fig. |4]shows the variation with temperature and den- 
sity of the number of correlated molecules N^J, in 
DBP, as obtained from Eq. ([TTJ in conjunction with Eqs. 
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12). N*IETa is found to increase when the tem- 



perature decreases or when density increases, as antici- 
pated. The number of dynamically correlated molecules 
obtained through this method is furthermore always 
slightly larger than those calculated through the other, 
less efficient, bound in Eq. d4b or ( 10 1. 



V. N corrA ALONG THE GLASS TRANSITION LINE AND 
OTHER ISOCHRONES 

From the 3-dimensional plot displayed in Fig. 4, one 
can follow the evolution of the number of dynamically 
correlated molecules along any chosen path in the (T, p) 
plane. An interesting choice of paths is provided by the 
isochronic, i.e. equal relaxation time, lines. The most 
commonly considered among such lines is the glass 
transition line, where TV, is defined by a given value of 
the relaxation time, say T a = 100 sec. It was reported 
in Ref. [16 1 for 4 glass-forming liquids and polymers 
that the number of dynamically correlated molecules 
is uniquely determined by the relaxation time. If this 
were true in general, N^JJa should be constant along 
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FIG. 4. 3D-plot for the number of dynamically correlated 



ucaiiyc 

molecules of DBP, estimated via Eqs. |6j, jllj and |12|, 

for a wide range of T and p. Nili.^ increases with increas- 
ing density and with decreasing temperature. The changes 
of color schematically indicate isochronic (constant relaxation 
time) conditions. 



isochronic lines. As shown in Fig. [5] (a), this is however 
not valid in the case of DBP. We find in particular that 
^corr 4 mcreases with pressure (and temperature) along 
the Tg line, at least up to the highest pressure available 
(the evolution appears to saturate or reach a maximum 
for the highest pressures or temperatures considered). 
The same systematic trend is observed for all isochrones, 
although the increase in becomes smaller as the 

relaxation time decreases [see Fig. |5](b)]. We stress that 
the pressure range covered in the present study on DBP 
is significantly larger than in Ref. [16]: here, we con- 
sider P from the atmospheric value up to 1 GPa (with a 
change in density of almost 25 %) whereas P varies from 
the atmospheric value to a few tens to a few hundreds 
MPa in Iflfel . 

Two remarks are worth making in connection with the 
above result. First, the stretching parameter f> of the re- 
laxation functions was shown to be constant along any 
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FIG. 5. -a Evolution of the number of dynamically correlated 
molecules along the glass transition line for DBP (Tg is 

defined for a dielectric relaxation of 100 sec), -b Same as (a) 
but for various isochronic lines (the value of the log 10 of the 
dielectric relaxation time in sec is given in the caption) 

isochronic line, and in particular along the T„ line, in 
several glass-forming liquids |25 28 1. The same obser- 
vation applies here to liquid DBP. Secondly, the isochoric 

fragility nip, defined as nip = ( ^ff^Npj^ ) * s a ^ so f oun d 

to be constant along any isochronic line, as a result of 
the density-temperature scaling of the relaxation time 
ET]|37| . The important consequence is that neither the 
stretching of the relaxation nor the dynamic fragility are 
correlated with the increase in the number of dynami- 
cally correlated molecules observed along the glass tran- 
sition line and other isochrones. 

As the stretching parameter j6 and the isochoric 
fragility m p are essentially constant, whereas Cy varies 
by 10 % along the glass transition line, it follows from 
Eqs. |llj | and (|sj that the variation of N^f^ 4 with pres- 
sure should be mostly due to that of the parameter R. 
This is indeed what is observed, as illustrated in Fig. [6] 
1 + R increases by a factor of almost 3 (and R by a factor 
of almost 4), which roughly corresponds to the increase 



FIG. 6. Variation along the Tg line of 1 + R, where the param- 
eter R quantifies the relative importance of density-induced 
versus temperature-induced contributions to N^L;,. As for 
the number of dynamically correlated molecules, 1 + R ap- 
pears to saturate or reach a maximum at the highest accessible 
pressures (and therefore temperatures) 



cut i 

The physical significance of this increase in R is that 
the fluctuation effects triggered by density become in- 
creasingly dominant over those associated with temper- 
ature as pressure and temperature increase. We stress 
that this should not be confused with the relative im- 
portance of density effects (versus temperature effects) 
in the slowing down of the relaxation. While R concerns 
the fluctuations around the average dynamics, the latter 
deals with the average dynamics (see the Introduction) 
and, in the scaling description of the relaxation time, it 
is measured by the factor Ta.p(T,p)x(p) [37]. For DBP 
along the Tg(P) line, Tap is equal to 0.12 at atmospheric 
pressure and decreases by about 25 % at high pressure, 
whereas, as already quoted, x varies from 2.5 to 4. As 
a result, the factor Tapx changes at the glass transi- 
tion from 0.30 to 0.36 over the same pressure range. At 
Tg(P), the influence of density on the average dynamics 
at constant pressure is thus significantly less than that of 
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temperature and increases by only 20 % under pressure, 
a result very different from that found for R in connec- 
tion with the dynamical heterogeneities. 

As the thermodynamic factor pTKjCy varies weakly 
with pressure along the glass transition line, the increase 
in R can be assigned to the change in the activation en- 
ergy parameter x(p) [compare with Eq. jl2|]. The lat- 
ter indeed increases from 2.5 to 4 along the Tg(P) line, 
which, once squared, explains the growth of R and of 
^corri- T ne P resent study therefore highlights the role 
of the deviation from simple power-law behavior of the 
bare activation energy e(p) in the variation of the num- 
ber of dynamically correlated molecules along the glass 
transition line. As a check, we have repeated for an- 
other molecular glass-forming liquid, ortho-terphenyl, 
the very same procedure as done here. This leads 
to a value of the number of dynamically correlated 
molecules that increases along the glass transition line 
with increasing pressure up to 400 MP a, but only by less 
than 20 %. Accordingly, it was found that the parameter 
x(p) can be taken as essentially constant over the same 
domain of pressure ETlEHl . 

Finally, we note that our finding (for DBP up to 1 GPa) 
of an increase of the number of dynamically correlated 
molecules, and as a consequence of the associated cor- 
relation length, along the Tg(P) line is inconsistent with 
the theoretical predictions made on the basis of the Ran- 
dom First-Order Transition (RFOT) theory ||26l |27| . In 
the RFOT theory, the dynamic length at T ? is predicted 
to be £ = 5.8 a 11261 |27| . where a is the diameter of an 
elementary component of the molecule or the polymer, 
referred to as a "bead". This prediction, which only de- 
pends on the relative value of the relaxation time com- 
pared to a microscopic time, is claimed to be valid for 
all glass-formers under any thermodynamic conditions. 
Some uncertainty do exist as to how beads should be 
defined and operational procedures have been devised 
to compare different systems 11261 |27| . However, when 
looking at the same liquid under different pressures, the 
bead should keep unchanged, which sidesteps the fuzzi- 
ness of its definition. The fact that N C orr,4 increases along 
the glass transition line then directly contradicts the 
claim of the RFOT theory (except if one thinks that the 



relation between length scale and number of molecules 
should change with pressure, an a priori groundless as- 
sumption). One may of course take the variation by a 
factor of slightly less than 3 which is observed observed 
for DBP as "within the noise" of the theoretical predic- 
tions and the experimental procedures. Nevertheless, 
considering the care with which the various sources of 
uncertainty have been estimated in our calculation, it 
rather seems as a fact to be seriously taken into account 
in further theoretical and experimental investigations. 
VI. CONCLUSION 



In this work, we have characterized the spatial extent 
of the dynamical heterogeneities in a glass-forming 
liquid through an experimental determination of the 
number of dynamically correlated molecules obtained 
from dynamical susceptibilities. We have focused on 
the fragile liquid dibutylphatalate (DBP), for which, 
from our own measurements and from existing ones, 
a very large domain of temperature and pressure is 
experimentally covered for both dielectric spectroscopy 
and thermodynamics. We have calculated different 
estimates of the number of dynamically correlated 
molecules N^J,, including a crisp one making use of 
the separate contributions due to fluctuations associated 
with density and with temperature. We have shown 
that the range of the spatial correlations in the dynamics 
of DBP grow when approaching the glass transition 
along any thermodynamic path. In addition, we have 
found that N^ r ^ varies along the glass transition line 
(by a factor of about 3), and more generally along any 
isochronic line, as pressure (or equivalently tempera- 
ture) increases, a result which is at odds with recent 
reports and theoretical predictions. The spatial extent 
of the correlations in the dynamics is further found to 
be uncorrelated with both the (isochoric) fragility of the 
system and the stretching of the relaxation function. 
This contradicts the naive view of the dynamics of 
glass-forming liquids that ties cooperativity fragility, 
stretching and growing heterogeneity altogether. 
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